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The quadrupole formula in four-dimensional Einstein gravity is a useful tool to describe gravita- 
tional wave radiation. We derive the quadrupole formula for the Kaluza-Klein (KK) modes in the 
Randall- Sundrum braneworld model. The quadrupole formula provides transparent representation 
of the exterior weak gravitational field induced by localized sources. We find that a general isolated 
dynamical source gives rise to the 1/r 2 correction to the leading 1/r gravitational field. We apply 
the formula to an evaluation of the effective energy carried by the KK modes from the viewpoint 
of an observer on the brane. Contrary to the ordinary gravitational waves (zero mode), the flux 
■ of the induced KK modes by the non-spherical part of the quadrupole moment vanishes at infinity 

and only the spherical part contributes to the flux. Since the effect of the KK modes appears in 
the linear order of the metric perturbations, the effective energy flux observed on the brane is not 
always positive, but can become negative depending on the motion of the localized sources. 
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INTRODUCTION 



Strong gravitational phenomena, such as gravitational collapse and black hole formation, have long been intriguing 
t-H \ subjects of General Relativity in the context of astrophysics. Likewise, studying strongly gravitating dynamics in 
frameworks of recent models of the universe in which we live on a 4D membrane (braneworlds) will promisingly open 
up new aspects of the models themselves [lj, 0, HI • One of the sites is naturally gravitational collapse occurring 
typically in the last stage of stellar evolution or during the early era of the universe (the latter corresponds primordial 
black holes), which would experience sufficiently high density and energy so that effects of extra dimensions might 
not be negligible. 

The model proposed by Randall and Sundrum (RS) 0,0] is a simple, but successful model that realizes a braneworld. 
i^l So far no serious drawback to the model has been known. (See, e.g., for a comprehensive review). In the model, 

£h . perturbative analysis of weak gravity reveals that the KK modes give the correction (oc r~ 3 ) to the Newton potential 
(oc r -1 ) for static isolated sources [a, ISM; Furthermore it is known that the gravity induced by relativistic stars 

^ . and black holes localized on the brane [ill Il2t Il3j | is well described by the lower dimensional 4D Einstein gravity on 
the brane. The induced gravity in the latter case is approximated by the black strin g: wh en the mass is sufficiently 
large, although the subject of the black hole in the braneworld is a little controversial |l4l fl5l| (see also, e.g. 0]). 

Our knowledge of strongly gravitating dynamical processes in the braneworld is, however, very limited. Even for 
a simple spherically symmetric collapse, the external gravitational field including Kaluza-Klein (KK) modes excited 
during the process has not been fully investigated. In 4D General Relativity, as a result of Birkhoff's theorem, 
all spherically symmetric vacuum spacetimes become the Schwarzschild solution even implying that any spherical 
collapse cannot radiate gravitational waves; only non-spherical dynamical processes produce gravitational waves. 
Their amplitude and energy loss rate are well described by the so-called quadrupole formula, whose prediction for 
the energy loss rate in the binary pulsars, e.g. PSR1913+16 and PSR1534+12, precisely agree with the observed 
orbital decay rate 01 ■ Thus the quadrupole formula is of interest to test alternative theories of gravity 0] . Even in 
higher-dimensional spacetimes, as far as we consider spherically symmetric collapse, the situation does not change 1191 . 
The higher-dimensional quadrupole formula in the asymptotically flat spacetimes has been recently discussed in [2(j 
(for even Z?-dimensions) . However, in realistic braneworld models such as RS, matter fields are confined on the brane 
and seemingly spherical collapse of such fields is not really so from the viewpoint of the bulk so that the Kaluza- 
Klein modes would be inevitably induced on the brane and the non-local effect from the bulk gravity complicates the 
problem. In considering homogeneous spherical dust collapse on the brane it was shown that outside the dust 
sphere the induced metric could not remain static and vacuum. Hence it is necessary and important to investigate 
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the dynamical process including dynamically generated KK gravitons, it will allow insight into the BH formation and 
astrophysical implications of the braneworld I2ll.l22ll23l .l24j. (See e.g. [25| for gravity wave perturbations using the 
black string.) 

The purpose of this paper is to have insight into physically interesting dynamical processes in the braneworld. 
Toward this direction, we investigate a gravitational field induced by an isolated dynamical system on the brane. We 
employ the linear perturbation theory of the RS model [M S to accomplish a self-contained analysis. We first derive 
a quadrupole formula which describes the (KK) gravity outside matter distribution by means of quadrupole moments. 
By applying this formula, we can successfully describe weak gravity far away from a source once quadrupole moments 
of an isolated system are determined. Based on the quadrupole representation of the induced metric, we also discuss 
effective energy that is transferred by the KK modes. This kind of energy flux corresponds in the linear order to the 
flux of so-called dark radiation in the braneworld cosmology [2(|. We find that the KK modes contribute to the flux 
only through the projected bulk Weyl tensor as an effective matter, but not through energy-momentum tensor of the 
gravitational field itself. Moreover, this effective energy density is not always positive, and in fact it becomes negative 
in the case of the homogeneous dust collapse. 

The paper is organized as follows. In Sec.[n] we review the linear perturbation theory of the RS model discussing the 
Green's function for the dynamical source, and derive the quadrupole formula for dynamical systems. In Sec. IIHI we 
apply the formulation to a spherically symmetric system on the brane. As a concrete example, we discuss spherical dust 
collapse and estimate the change rate of the effective energy associated with the KK modes. General non-spherically 
symmetric cases are also discussed in this section. 



II. LINEAR PERTURBATION IN THE BRANEWORLD 



A. Green's function 



Let us begin with metric perturbation in the RS model in order to discuss weak gravitational field generated by an 
isolated dynamical system on the brane. We follow the procedure given in (Notation in this paper is the same as 
in §.) Denoting metric perturbation in the RS gauge as h^, the 5D metric is given by 

ds 2 = — [dz 2 + r } „vdx fi dx v ] + h^dx^dx", (1) 

where 7y M „ is the 4D Minkowski metric and I is the curvature radius of anti-de Sitter (AdS). In this gauge linearized 
5D Einstein's equation becomes 



dz 2 z dz z 2 



h ltv = -2KEv tt ,-5(z-£), (2) 
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where the source term is given by 

( 1 \ £ 
%iiv = [Tpv- -Tp^ T ) + 2k C 7/w = V- (3) 

Here is the energy-momentum tensor of matter on the brane at z = I and £ 5 expresses location of the brane in the 
RS gauge, which is related to the matter source on the brane @. 4D Newton's constant is G4 and we use re = 
Using retarded Green's function the formal solution of Eq. |J3J is obtained as 

h flu = —2K / d 4 x'G(x, z,x',£)T,^ u (x'), (4) 



where G(x, z,x' , z') is decomposed into two parts; 

G(x, z, x , z ) = G zcro + Gkk- (5) 

G zcro is the usual Green's function of massless scalar fields and Gkk is that of massive Kaluza-Klein modes. Cor- 
respondingly perturbed metric is written as h^ v = h z ^° + h^, where we note that for the RS single-brane model 
the zero mode corresponds to the linear perturbation of 4D Einstein gravity @ . The Fourier components of the KK 
modes are given by 



DC: 



d 3 k u m {z)u m (£) o4k . (x _ x ') 



Gkk( W ,xV) = - / dm/ (6) 
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where the mode function is u m (z) = \ mL J ^ mt )J 2 ^ mz ) Y 1 {mt).j- 2 {mz) ^ j n ^ g p a p er jjat d eno t es Fourier components. 

We are interested in the case that matter is non-relativistic and its density is low, hence the induced gravity is well 
described by linear perturbations. Then only low-frequency modes to <C 1/i are relevant. (The mi 3> 1 case was 
investigated in jUj.) Expanding the mode function about ml <§; 1, one gets 

~ f°° mP f f\ 2 f H 3 k pik-(x-x') 



u 



1 f°° , ml I 



inR J 

i 3 ( iuj i 



2 



dm^(-) e iR ^^ + 0(m 3 £% (7) 



e iu>R 



8nz 2 \R 2 R 3 , 

where R = |x — x'|. Here we are interested in the gravitational field on the brane so that we have set z ~ I. Note 
that very massive modes with m > to are exponentially suppressed by the factor e -RVm 2 -uj 2 j n ^ e re gj on f ar f r0 m 
the source, while only the KK modes with m < lu contribute to the m integral. In the static case (lu = 0) the first 
term vanishes and the second gives the well-known KK correction to the Newtonian potential, which is proportional 
to l 2 /R 3 [1 S @. In the dynamical case (uj 7^ 0), however, the first term survives. This term corresponds to 
dynamical KK modes. This kind of difference of falloff between dynamical and static terms is known to occur in 
higher dimensional Green's functions [2(j , whereas for the 4D massless Green's function both the static potential and 
the wave form (dynamical term) fall as r _1 toward infinity. 

Compared with the Einstein gravity produced by G ZGr o, which is proportional to l/R, the dynamical KK field 
seems to be suppressed in the region far from the isolated matter. However, the radial dependence of the metric for 
the KK modes is not necessarily suppressed compared to that of the Einstein gravity, as we will see explicitly below 
employing the areal radius for the spherically symmetric case. Thus in this sense, the dynamical KK modes could 
produce "leading order" contributions, and then it is interesting to study various effects induced by the dynamical 
KK corrections. 

Before closing this subsection, it is worthwhile to note that the above dynamical KK corrections could be interpreted 
into the context of conformal field theory via the AdS/CFT correspondence in the braneworld. The correspondence 
has been ever discussed to investigate for example following aspects of the braneworld gravity; 1 / R 3 correction to the 
Newtonian gravitational potential, homogeneous cosmology, and tensor-type perturbations on a FLRW brane 0, 
l2Sl EflllsnL l3ll| . It is quite interesting to consider if this correspondence also persists in dynamical cases. 

For dynamical sources (w ^ 0), the CFT one- loop correction ^(p) to the graviton propagator |22| are reduced in 
the Fourier space to 



(2 r oo 



— oo 



e i 

167T V r 2 r 3 



&ke lkr k{\n{k -u-ie) + In (A + lu + ie)} (8) 



where ^(p) ~ ^lnp 2 and p 2 = —(lu + ie) 2 + k 2 . This means that in dynamical cases the KK Green's function of 
Eq. JJJ can be rewritten by the CFT correction ^(p 2 ). Therefore it is possible that the dynamical external field 
discussing below could be interpreted via CFT description. 



B. Quadrupole formula for KK fields 



In this subsection we shall write the metric perturbation of the KK excitations in terms of quadrupole moment 
of matter distribution. The quadrupole formula is very common in General Relativity. In 4D General Relativity 
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corresponding to the zero mode, the metric perturbation is given under an appropriate gauge 32] by 



+0 (r- 3 ), 



(9) 



W»(t,x)=^i 4 



r dt 2 10 



I^(t-r)+0(r-<) 



where lJ } T is transverse-traceless components of quadrupole moment and J y = J d 3 x(x i T^° — x^T M ) is total angular 
momentum. Thus in these leading terms only the transverse-traceless components are time-dependent, which describe 
gravitational waves. The other terms describe total mass M and total angular momentum Jy , which are conserved 
quantities (time-independent). 

Let us now return to the KK modes. From Eq. the metric perturbation of the KK excitations on the brane at 
z = £ is 



V ( W) x)~-2G 4 ^ 



dVe^ 



ILU 1 



T»u ~ ^V^T ) (w,x'). 



(10) 



We note that h^ v is obtained by the gauge transformation to the Gaussian normal coordinates, and the gauge 
transformation ^ in the source term of Eq. does not contribute to the KK modes Q . 

To obtain the quadrupole representation of the metric perturbation, we rewrite the integral of the energy- momentum 
tensor in terms of quadrupole moment. By energy conservation, 9^T A " y = 0, we have the following relation j3^|, 



J d 3 xf « = J d 3 x[d fc (f fc V) - x l d k f k3 \ 



d 3 xx 4 T 0j 



iui 

"T 

uj 2 
' 2 



d 3 x[a fe (f ofe xV) - x l x>d k f ' 
, ,2 



U/vl 



(11) 



d 3 xxVT 00 



where P J is quadrupole moment and k are spatial indexes. At this point, we assume w|x'| <C 1 in the slow-motion 
approximation that the internal velocities of sources are small and thus the typical source radius is much smaller than 
the wavelength l/u>. At large distances from the source, i.e., r = |x| 3> x', it is possible to make multipole expansions 
up to the order of quadrupole moment: 



d 3 / JuR \^__J_\ T t3 

R 2 R 3 



_ — 1 foj 

R 2 R 3 

^L-—\ T" 
R 2 R 3 ' 



(— iuj) 3 (— iuj) 2 



\3 i s,_,\2 



I l 3{iu)e l 



(-«") , , N + 3 (-w) 



^2^(c) - i 

r J 2 
(— iw) 3 (— io;) 



(-iuj) 3 , 2 (~^) 2 



r 



3r 2 



3r 2 



.(-iw) 



5 1 



(12) 



Here in the rest frame / d 3 xT 00 = M and J d 3 xT 0i = 0, where M is the total mass of the system. Substituting these 
into Eq. I|10fl . we obtain the metric perturbation expressed by the quadrupole moment, the results in the real space 
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Gd 



r 2 ch 3 



where / = / 



^ K (t,x) 



and / , , = /, , 



Gd 2 



3r 3 



Gd 



J%2 



r 2 dt 3 



_3_ d^ ^_ d_ 

r 3 dt 2 r 4 d£ J r 

J_ _1_ d^ 

3r 2 di 3 + r 3 di 2 

d 2 " 



— Iij{t 



(13) 



2 d 
r 4 dt 



2 



J(i-r)*« 



1 



dt 2 



\I8ij is the trace-free part whose components will vanish if the source is spherically 



symmetric. The terms of 0(l/r ) which are independent of the quadrupole moments are static KK corrections. If 
the matter source is dynamical, then the dynamical terms of order 1/r 2 arise. 

The result obtained here provides a formula that describes the far exterior gravitational field induced by general 
matter sources localized on the brane. It shows explicitly that in the braneworld even spherical sources can make the 
exterior gravitational field dynamical. (See Q for the case of homogeneous dust collapse.) 

Although the linear perturbation is used to represent weak field approximation, it does not imply the lack of 
possibility to describe strong-gravity systems (e.g., a self-gravitational system). Because for asymptotically flat 
spacetimes the gravitational field will become so weak to be dealt as linear perturbations at sufficiently far region, the 
linearized Einstein e quat ion can be thought of as an "exact" equation if other nonlinear terms are taken as a source 
term of the equation |32l l34| . In other words, even if the interior is strongly gravitating system which needs nonlinear 
treatment, its far exterior field can be described in the same manner as for a weakly gravitating system except the 
definition of total mass or quadrupole moment. (See, e.g., |35).) 



III. APPLICATIONS 



A. Spherically symmetric case 



Now we consider a spherically symmetric source, of which consequences are non-trivial in the braneworld. By virtue 
of spherical symmetry, the quadrupole moment becomes 



1 



lijit) = M(t), I{t) = I k k = / d 3 x r 2 p 



(14) 



KK 



where the trace-free part -t ^ vanishes. In the spherical case the spatial components of the metric perturbation h 
are proportional to &y so that the metric can be expressed in the isotropic coordinates. The expression becomes 
concise after being transformed into the area coordinates by following gauge transformation 
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Gd 2 
18 



7 (3) + » /(2) + » 7 (i) 



MGd 2 Gd 2 



3r 2 



18 



17(3) + ±7(2) + ±7(D + ±7 



d" 



where 1^ = ~ r )- The non- vanishing components in the area coordinates are 



^ K (*,x) 



AMGd 2 Gd 2 



3r 3 

KK/, . 



2MGd 2 Gd 2 



1/(4) + jl/(3) _ i_j(2) _ IH7W _ 1^7 

I/( 4 ) + ±7(3) + ^7(2) + 30j(i) + —I 



(15) 



(16) 



Compared with 4D General Relativity @ , a new term proportional to 1 jr appears in the radial component /iJ5 r K in 
this coordinates which claims that the KK modes may contribute to energy evaluated at infinity. 
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As mentioned above, the expression of the metric perturbation depends on the gauge. We show the contribution 
to the energy using the Einstein tensor which is a gauge-invariant quantity in this order. The effective 4D Einstein 
equation on a vacuum brane is given by 



(17) 



where comes from 5D Weyl tensor. Since now we have the induced metric on the brane, we can directory obtain 
the effects from the bulk. Straightforward calculation of LHS in the linear order leads to 
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Gd 2 



2MGd 2 
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ZMGd 2 
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3^ 



20 
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10 

9^ 



j(3) 
j(3) 



20 

10 

3^ 



j(5) 



9r 3 



j(4) 



3r 4 



.j(3) 



r(2) 



10 

37s 



_25 

3^6 



(i) 



j(2) 

J(l) 

25 
3^7 



40 
3^ 

10 

3^ 



40 
3^ 



(18) 



/' 2 



10 

3^6 
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It is explicitly shown that £^ satisfies the local conservation law in the linear order, <9 M £ MI/ = 0, and the trace-free 
condition, = 0. Since can be regarded as an effective energy- momentum tensor on the brane, we can discuss 
the effective energy flux of the KK modes corresponding to £ A ,„. The total energy flux of the KK modes, £kk> is 
evaluated at future null infinity on the brane (r — * oo, t — r = const.) as 



dE 



KK 



1 



dt 



8irG 4 dt 



(-£ tt )d^ = --±- 

y 07T(_r4 



18 



:/ (5) 



(19) 



where E KK is the KK part of the energy of the system. Thus this is suppressed by I 2 . We should notice that Eq. Ijl9(l 
is not quadratic in the quadrupole moment as contrasted with the well-known expression for the 4D gravitational 
radiation. This is due to the fact that £ M „ exists in the linear order of metric perturbation. The sign of dE KK /dt is 
indefinite, in other words, the KK modes do not always carry positive energy away in contrast to the zero mode |40| . 
We will show an example in the next subsection. 

While been considering £^ v as the energy of the KK modes, in usual discussions of the energy of gravitational waves 
we use the energy-momentum pseudo-tensor of gravitational fields, which is quadratic in the metric perturbation |32l 
l3fij . After a straightforward calculation of the Einstein tensor to the second order, ^G^, we find 



{2) G, V ~ O 



(20) 



Thus the contributions are higher order, meaning that this kind of energy for the KK modes cannot reach infinity. 
Consequently the KK modes carry its energy to null infinity on the brane not through the quadratic order (gravitational 
energy) but the linear order of like an effective matter flux. We will observe another aspect of this nature in 
Appendix, in which we relax the assumption of the slow motion and it turns out that £^ v has a flux part explicitly. 



B. Spherical dust collapse 

As a concrete, simple example of dynamical systems we shall consider homogeneous dust collapse under spheri- 
cal symmetry. This model is well known as the Oppenheimer-Snyder model in the 4D case. As our attention is 
concentrated on the weak gravitational field induced by non-relativistic matter, we deal with the matter source as 
Newtonian. We denote the mass, radius, and density of the dust cloud by M, R, and p, respectively. The equation 
of motion of the dust sphere is given by 

dR\ 2 r s , . 3M 

&) *i> ^ = ^rWy (21) 

where rs = 2G 4 M is the 4D Schwarzschild radius. The quadrupole moment of this system is I(t) = ^MR 2 (t). Here 
we have neglected the KK corrections to the equation of motion and the time derivative of quadrupole moment, 
because they only give higher order corrections in the resulting metric h^^. 
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Substituting these into Eqs. I|16fl and (|19fl . we obtain 

I I /r s \3 



h KK 
and 



Lkk - ~6G7r5 (?) ' ' (23) 

These results indicate that the exterior metric of the collapsing dust is manifestly non-static and the total flux of the 
KK modes is negative. This means that the total effective energy of the system for an observer at infinity on the 
brane increases during the collapse. This is the direct consequence of the negativity of the effective energy density 
of E^u- The gravitational fields are expressed using the retarded Green's function, so outside the sources any flux is 
outgoing, in other words, depends on retarded time t — r. If the energy density of the flux is negative, the total energy 
on the brane will increase. This kind of negative energy density for £ ^ v is general in the braneworld. For example, 
the exact solution of the localized black hole in the 4D RS model [^3 or the KK correction to the Newton potential 
for a static source 8] give the negative energy density. Also in the brane cosmologies the KK modes behave as dust 
with a negative energy density on the brane |38| . 



C. Non-spherically symmetric case 



We shall now briefly discuss non-spherical cases. The terms including -f jj remain in Eq. 1)13(1 and give additional 
terms for £ M „ other than those in Eq. I(18|l . The relevant component to such energy flux is £t r which is evaluated as 



rnon 
c tr 



G/xV (5) 3 

3r 2 r 2 * V >> 



(24) 



This leads to the energy change rate due to the non-spherical terms. However, it turns out that the rate is zero, 



dE° 



dt 



-t t0 an 



3 °v 



0, 



(25) 



and that the non- spherical part of the KK modes cannot carry the energy. This means that for general dynamical 
processes such as particles infalling into a black hole or binary stars, only the spherical part is relevant to the total 
flux of the KK modes. This is remarkable contrast to the emission mechanism of ordinary gravitational waves (zero 
modes). 

We shall compare the energy flux of the KK modes with that of the zero mode. The total energy flux of the zero 
mode is given by 



T _ ^i_t(3) .(3) 

-^zcro — r TiA ~r ; 
5 



«7 



and we estimate 



L 



KK 



£ £ (cT) 
G 4 MR R 



10 



-16 



/ f/0.1mm 

\G 4 M/M 



3/2 / l/0.1mm \ 1/2 
V R/lOkm J 



(26) 



(27) 



where M,R and T are characteristic mass, source radius and time-scale, respectively, taken to be I <C R < cT ~ 
y/ R 3 /G 4 M . Thus the relative amount of energy carried by the KK modes is negligible in general astrophysical 
situations. 



IV. SUMMARY AND DISCUSSION 



We have discussed the gravitational fields surrounding isolated dynamical systems in the RS model and obtained 
the quadrupole formula for the KK modes. We found that in dynamical cases the leading term of the KK correction 
for the Green's function is of order 1/r 2 rather than 1/r 3 , which is observed in the static case. By applying the 
quadrupole formula, we have considered the effective energy flux of the KK modes. Contrary to the fact that ordinary 
gravitational waves (the zero mode) are generated by non-spherical parts of matter distributions, the effective energy 
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flux of the KK modes arises only through the spherical part. This implies that it is difficult to generate the flux of the 
KK modes by a circular motion such as binary pulsars. Moreover, even for spherical collapses which cannot radiate 
the zero mode, the radiated KK modes are suppressed by £ 2 as well as static cases @, El • 

For spherical dust collapse we have explicitly derived the exterior metric on the brane and evaluated the effective 
energy flux during the collapse. In our treatment, the KK modes contribute to the total energy of systems on the 
brane through the projected bulk Weyl tensor, but not through the energy-momentum tensor of gravitational 
field itself. This result is a consequence of the fact that in the 4D picture of the braneworld the KK modes induce 
the energy-momentum as an effective matter in the form of —S^v 

In the linear order, satisfies the local conservation law as the energy-momentum tensor of matter fields satisfies, 
and hence it has independent conserved charges at this order. In this sense the effective energy flux of the KK modes 
is irrelevant to the energy loss via gravitational wave radiation. To derive the correction to the energy loss rate due 
to KK gravitons, we will need to go to second-order perturbations |4l|(see also, e.g., f° r ordinary gravitational 
waves). 

In nonlinear regime the KK modes couple to the matter distribution due to the relation V M (£ M „ — x ^nv) = 
where 7r M „ consists of the energy- momentum tensor T^ v |3flj . and it does not have a separate conserved charge. Thus 
the next-order perturbations will affect the energy of matter, and it will reveal us how the ordinary gravitational 
wave forms and the energy loss rate are affected by the KK modes. It is interesting to derive a quadrupole formula 
as a result of the second-order interaction and we expect to provide new tests for the braneworld. The analysis of the 
nonlinear interaction and energy loss rate will be discussed in a forthcoming paper |43j . (A brief discussion is given 
inRef. 0) 

The flux of the KK modes works as an effective matter flux, namely a flux of the dark radiation. This type of 
energy flux is not necessarily to be positive. In fact, as we have observed in the case of the dust collapse, the KK 
modes have negative energy density. As a result the total charge associated with £^ v increases after the collapse. 
This is related to the fact that the 5D Weyl tensor in the effective Einstein equations appears in the linear order 
of perturbations so that the sign of £^ v depends on that of the metric perturbation. 

We expect that the results can be understood in terms of the AdS / CFT correspondence because it was confirmed 
that the KK Green's function for dynamical processes corresponds to the CFT correction. It is interesting to give 
physical meaning to the negative energy flow via the CFT picture [44| . 
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APPENDIX A: BEYOND SLOW-MOTION APPROXIMATION 

In this appendix we relax the assumption of the slowly moving matter sources used throughout t his p aper. Then 
one can argue the gravitational wave from, for example, head-on collision of relativistic particles [2(], l36j . 

Let us first recall our treatment within the slow-motion approximation (o->|x'| <g; 1). We could write the integral in 
Eq. 112|) in term of the quadrupole moments as 

J d 3 x'e iw|x - x,| T ij (w,x') ~ e lujr J d 3 x'T y '(tj, x') = -^-rJ(w)e <u "", (Al) 

where |x'| <§C |x| = r and we set e~H x > ~ 1. However, in more general cases, we cannot neglect the phase factor 
e -iu\x \^ an( j foe above integral is estimated as 

f d 3 x'e iw|x - x,| T^(cj,x') ~e wr f d 3 x'e~ M T^T^(w, x') 

J \ (A2) 

=e lujr / d 3 x'e- 4k ' x 'T^(w,x') = T^(w,k)e Mr , 



where k = wx/|x| and the last equality is the Fourier transformation with respect to spatial coordinates. Note that 
the conservation law yields k^T^" — where k^ = (u>, k) and k^k^ = 0. 
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From Eq. I|10|l . the leading order of the KK modes without the slow- motion approximation becomes 

h™ =- 2 *J d 3 x'G KK (t^ - ^V r ) (^' x ') 

/ , \ (A3) 
~ - 2G 4 £ 2 ^e^ r (f MV - - V ^f J (a;, k) + ©(r" 3 ), 



and the corresponding projected Weyl tensor is 



^ = ^!^Tfc^ + 0(r- 3 ). (A4) 



This explicitly indicates that the leading term of £ M „ behaves as a null fluid. 
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